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Abstract 

We show that the time evolution of the wave function of a quantum 
mechanical many particle system can be implemented very efficiently on 
a quantum computer. The computational cost of such a simulation is 
comparable to the cost of a conventional simulation of the corresponding 
classical system. We then sketch how results of interest, like the energy 
spectrum of a system, can be obtained. We also indicate that ultimately 
the simulation of quantum field theory might be possible on large quantum 
computers. 

We want to demonstrate that in principle various interesting things 
can be done. Actual applications will have to be worked out in detail also 
depending on what kind of quantum computer may be available one day... 

1 Quantum Computers (QCs) 

Quantum computers are still imaginary devices, but it is hoped that even- 
tually the technical problems involved in their realization can be overcome 
[0, [|, ||, ||, ||, ||, @ ■ Quantum computers could solve some problems much faster 
than conventional computers. Most prominently, Peter Shor (1994) has given 
a "quantum algorithm" for factoring large integers in polynomial time f8|(see 
also §). 

A /-bit quantum computer may be thought of consisting of I two-state sys- 
tems. Computations would be carried out by inducing unitary transformations 
of a few at a time of these quantum-bits (qubits). This may be done by exterior 
fields controlled from the outside. Upon observation at the end of the "unitary 
computation" the quantum computer would collapse into a state where each 
qubit is either a or a 1. The art of finding a "quantum algorithm" is to 
extract useful information from a few such runs. 

The main technical problem in realizing QCs is, to prevent unwanted in- 
teractions with the environment during the calculation. Such interactions can 
cause decoherence, such that the QC is no longer in a pure quantum state. Here 
we consider an idealized QC which doesn't suffer from such problems. 
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2 Simulating Quantum Systems 



General ideas about using specially designed quantum systems to simulate other 
quantum systems have been published, e.g. by Feynman ||. I present here an 
actual implementation of the simulation of quantum mechanical many particle 
systems on a general purpose QC q. For every degree of freedom of the system 
we need a Z-bit quantum register. This allows to store the whole (discretized) 
wave function as amplitudes of the "classical" states in the quantum computer. 
Already to store so much information would be impossible for a conventional 
computer. It is then shown that the time evolution for short time steps amounts 
essentially to a Fourier transformation which can be carried out very efficiently 
on a quantum computer J§|, [l^, |ll[ . 

2.1 Quantum Mechanical Particle in 1 Dimension 

First we discretize the wave function and impose periodic boundary conditions: 

a n = ip{n Ax) a n+N = a n . (1) 
After proper normalization we store these amplitudes in a i-bit quantum register: 

JV-l 

= ]T a n \n) N — 2 l . (2) 

Where \n) is the basis state corresponding to the binary representation of the 
number n. For short time steps At the Greens function of the Schrodinger 
equation is approximately: 

G( Xl ,x 2 , At) = k e*9 ^^-iV(x 1 )At (3) 

Applying this to the amplitudes a n is equivalent to transforming the basis states 
with the inverse transformation. Using proper normalization to make the trans- 
formation unitary, we get: 

1 N ~ 1 _ I 2 2 

In) -> -== T e-^ ( """At A ' A -> At In') . (4) 

, „ 

v n'—0 

The crucial observation is that this corresponds to two diagonal unitary matrices 
plus a discrete Fourier transformation: 



, N-l 

1 e _ if a^d+fl^,, Ax) At ^ ^m^l f --mi^ 

N > r 



jv — i 

^ e im^^- L- if 2^1^ , ( 5 ) 



dia S° nal Fourier trafo dia S onal 

'Independently and at practically the same time as me, Steven Wiesner Jl^ has published 



a preprint concerning this, using a slightly different approach 
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To be able to use the fast Fourier transformation algorithm (FFT), we chose 
the parameters such that: 

AX 2 „27T 

m— — = A — , A integer . (6) 

At N 6 w 

Provided the potential V has been made periodic (with period N), this also 
makes the whole phase factor periodic in n and n', as is necessary for consistency. 
It is rather straightforward to implement the FFT algorithm on a quantum 
computer ||, [Io| |ll[ such that it uses only about l 2 /2 computational steps (local 
unitary transformations) . The diagonal unitary transformations of the type 

|n) -> e icF ^ \n) , (7) 
can be done with the following succession of steps: 



|n,0)-|n,F(n)) ; \F(n)) - e* cF ("> \F{n)) ; \n, F(n)) - |n, 0) . 

(8) 

In the first and last step the vectors represent two quantum registers. The first 
step corresponds to parallel classical calculations. It can be carried out on a 
quantum computer with a somewhat higher cost than just one such calculation 
on a conventional computer (P. Shor ||). As all transformations are unitary 
they can just as well be inverted, so that the last step is also possible. The 
second step is really a transformation of the form 

\n)^e icn \n) , (9) 

which can be carried out in a straightforward manner in I steps by transforming 
each bit separately by an appropriate transformation. Note that of course all 
these operations are carried out in parallel on all basis states which is why this 
is often called " quantum parallelism" . 



2.2 Many Particles and Field Theory 

The generalization to many particles is straightforward. For n particles in 3 
dimensions we need 3n quantum registers. Besides the above steps we need 
transformations corresponding to the coupling between particles and between 
the different degrees of freedom of the same particle. They are diagonal unitary 
transformations acting on several registers, e.g.: 

\n, n', n") -> e lcF{n ' n ' \n, n', n") . (10) 

Such a transformation can be carried out analogously to eq. (j8j) . For identical 
particles the initial state of the quantum computer has to be chosen symmet- 
ric resp. antisymmetric. Also the generalization to more general Hamiltonian 
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operators than in the previous paragraph is possible (e.g. with coupling to a 
magnetic field). 

Also quantum field theory could be simulated, e.g. by discretizing the field as 
in lattice gauge theory (for a standard text see e.g. p5|). For a bosonic quantum 
field this would be equivalent to having 1 particle at each lattice point. At any 
rate, applications like lattice QCD would require large quantum computers with 
thousands of quantum bits, which are not going to be realizable soon. 

2.3 Fermionic Field Theories 

Fermionic quantum fields pose some problems, as their functional formulation 
(path integrals, wave functionals,...) involves anticommuting Grassmann num- 
bers. There is no wave functional in terms of the usual complex numbers, so 
we have to look for a representation of the fermionic field operator-algebra on 
some other Hilbert space. One possibility is the Fock space, where we have 
for the n-particle sector a totally antisymmetric wave function of n variables. 
Another possibility is to give the occupation number for each particle species at 
each lattice point. For fermions this can only be or 1 and thus requires only 
1 qubit. 

I expect, that further difficulties which may arise with quantum fields can 
in principle be overcome, possibly along the same lines as it is done in lattice 
QCD. 

3 Other manipulations 

3.1 Simulating a decay to obtain the ground state 

Often one is interested in the ground state of a quantum system, be this the 
ground state of a molecule, the vacuum of a field theory or a hadron stable under 
the strong interaction. I propose here to obtain the ground state essentially in 
the same way as it happens in reality, namely by letting some initial state decay. 
Decays happen because the instable particle is coupled to some other quantum 
system, like the electromagnetic field. It is important that this quantum system 
has a lot of different energy eigenvalues such that it can absorb the energy from 
various transitions in the decaying system. For the simulation of a quantum 
computer we may take a collection of 2-level systems with energy gaps AE = 
Eq 2~ n , n = 0. . .1 such that the whole system has practically a continuous 
energy spectrum. The coupling to the quantum system of interest may be done 
in various ways. 

During the simulation of the decay, the auxiliary "energy-drain" system 
must periodically be reset to its ground state such that it is able to absorb more 
energy. As the study of simple systems shows, this does not necessarily in every 
step, but only in the average, lead to a reduction of the energy expectation 
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value of the system of interest. For actual implementations one may have to 
think about an adequate "resetting" strategy, e.g. depending on the state the 
auxiliary system is found in upon observation just before each resetting. Also 
one may want to reduce the coupling towards the end of the decay simulation 
in order not to disturb the ground state of interest too much. 

Actually it may be enough to simulate the decay until the few lowest lying 
states (which are of int erest ) have a large amplitude. Then one can follow the 
procedure of paragraph 4.2 to obtain one of the energy eigenstates. 



3.2 Putting a wave function on a Quantum Computer 

Here I sketch a method which in principle allows us to store a given wave function 
as an initial state on a QC. I demonstrate this for the case of a scalar quantum 
mechanical particle in 1 dimension. The generalization to multi dimensional 
wave functions is then rather straightforward. Here I show only how the absolute 
value of the wave function can be stored in a quantum register, as above I have 
already shown how the complex phase can later be introduced. 

So we want to carry out the following transformation on a I bit quantum 
register which initially is in the state |0): 



|0> - fc £ 



where ip(x + L) = if)(x) 



(11) 



and k is an appropriate normalization constant. The idea now is to split the 
norm of the initial state (which is 1, of course) I times such that in the end 
we get the 2 l contributions for the single basis states (Imagine a binary 

tree with I levels and 2 l leaves.) For this to work out correctly each split has to 
yield the right ratio. For this we must have an algorithm which can compute 
the integrals 



I ik = 2 \f{x)\ 2 dx , for k = 0...2 i -l and i = l...l . (12) 
J XL 

Now the first split is realized by a real 0(2) rotation acting on the most signif- 
icant qubit of the quantum register, thus: 

|o,...)^yi^|o,...> + v^7|i,...> . (13) 

Such an 0(2) rotation of a qubit can be carried out analogously to the rephasing 
of an amplitude (see eq. ((?])) once the rotation angle ip has been calculated: 

2 Note that the norm of mutually orthogonal substates is a useful notion in quantum com- 
puting because it is conserved under unitary transformations. 
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sin ip — a/ Ji,i 



(14) 



The subsequent calculations of angles and rotations are then carried out in 
quantum parallelism, thus: 

\bi,bi-%, ...k-i,0, •■■) -> cos(vj) |6j, 6j_x, .. 0, ...)+sin(y>) \b h bi-i, ...h-i, 1, ■■•) 

. (15) 

with sin(( y 9) = ^//^j+x.fe and fc = + 2 + ... + 2 % b l . 

So this happens simultaneously for all values of the 6's (resp. of k). 

The method given here to input a wave function to a QC is rather ambitious, 
as essentially it evaluates the function at every (discretized) point x. In practice 
a simpler sceme may do. 

4 Obtaining Results of Interest 

4.1 Measuring positions of particles, resp. field strengths 

One can think of various ways how to obtain quantities of interest from a quan- 
tum computer simulation of a quantum system. Say we have some state on the 
QC which we want to analyze, like the vacuum or a 1-particle state of a field 
theory. One can now make a measurement of some field strength simply by ob- 
serving the corresponding quantum register. By repeatedly observing some field 
values in different points one may in several runs of the QC get a statistical esti- 
mate of n-point functions. In particular, one may obtain the correlation length. 
To get n-point functions for space-time points with a timelike separation, one 
has to continue the simulation of the time evolution between observations of 
single quantum registers. 

4.2 Measuring arbitrary observables 

To study the measurement problem in quantum theory, J. von Neumann has 
proposed an idealized simple interaction [m| between the quantum system 
and the part of the measuring apparatus that directly interacts with the system. 
For convenience this part of the apparatus can be chosen to be equivalent to a 
quantum mechanical particle in 1 dimension or, more precisely, its Hilbert space 
is spanned by the basis vectors 

\x) x real, and X\x) = ^l^) ■ (16) 

The time evolution during the measurement of an observable A shall then be 
given by 



G 



U(t)\V a )\x) = \V a )\x + k-a-t) , with = a|* ) , (17) 



where |W a ) is the state of the system. One easily verifies, that the Hamilton 
operator H = kPA leads to this time evolution. Here P is the usual momentum 
operator for the imaginary quantum mechanical particle, thus [X, P] = ih ■ 

Now we simply have to implement this time evolution on the QC. This can be 
expected to work out for any reasonable observable represented by the hermitian 
operator A which can be written in terms of the fundamental observables and 
their conjugate momenta. 

Then, after simulating this "1. stage von Neumann measuring process", we 
can observe the auxiliary quantum register and should find a number propor- 
tional to an eigenvalue of A. By repeating this, we get the spectrum of A and 
the weight of the various A-eigenstates in the original system state. Also after 
each observation the system will be in an eigenstate of A. Thus this is a way of 
obtaining e.g. the energy eigenstates. 



5 Conclusion 

I have demonstrated that on a large enough quantum computer various quantum 
theoretic quantities of interest could be calculated, which may be too hard to 
compute by a conventional computer. This is in particular true for strongly 
interacting field theories like QCD. However simulating quantum fields would 
require large quantum computers with many thousand qubits, which we don't 
know whether they can ever be built. Furthermore, here I have assumed to have 
an ideal QC which doesn't suffer from decoherence and can carry out unitary 
transformations with sufficient precision (as these are really analog operations). 
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